SARS-CoV-2 introductions to the island of Ireland (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had an unprecedented impact on the people of Ireland as waves of infection spread across the island during the global coronavirus disease 2019 (COVID-19) pandemic. Viral whole-genome sequencing (WGS) has provided insights into SARS-CoV-2 molecular mechanisms of pathogenicity and evolution and contributed to the development of anti-virals and vaccines. High levels of WGS have enabled effective SARS-CoV-2 genomic surveillance on the island of Ireland, leading to the generation of a sizeable data set with potential to provide additional insights into viral epidemiology. Because Ireland is an island, accurate documentation of travel rates to and from other regions, both by air and by sea, are available. Furthermore, the two distinct political jurisdictions on the island allow comparison of the impact of varying public health responses on viral dynamics, including SARS-CoV-2 introduction events. Using phylogenomic analysis incorporating sample collection date and location metadata, we identify multiple introduction and spreading events for all major viral lineages to the island of Ireland during the period studied (March 2020-June 2022). The majority of SARS-CoV-2 introductions originated from England, with frequent introductions from USA and northwestern Europe. The clusters of sequences predicted to derive from discrete introduction events ("introduction clusters") vary greatly in size, with some involving only one or two cases and others comprising thousands of samples. When introduction cluster samples are mapped sequentially by collection date, they appear predominantly in previously affected or adjacent areas. This mirroring of the phylogenetic relationships by the geospatiotemporal propagation of SARS-CoV-2 validates our analytic approach. By downsampling, we estimate the power to detect introductions to Ireland as a function of sequencing levels. Per capita normalisation of both sequencing levels and detected introductions accounts for biases due to differing sequencing efforts and total populations. This approach showed similar rates of introductions for all major lineages into Northern Ireland (NI) and Republic of Ireland (RoI) with the exception of Delta, which was higher in NI which is likely attributable to higher travel per capita. However, there were similar rates of Delta infection within NI and RoI, suggesting that although travel restrictions will reduce risk of introducing novel variants to the region, they may not substantially decrease total incidence. Our generalisable methodology to study introduction dynamics and optimal sequencing levels will assist public health authorities to select the most appropriate control measures and viral sequencing strategy.

MAJORA: Continuous integration supporting decentralised sequencing for SARS-CoV-2 genomic surveillance (preprint)

Genomic epidemiology has become an increasingly common tool for epidemic response. Recent technological advances have made it possible to sequence genomes rapidly enough to inform outbreak response, and cheaply enough to justify dense sampling of even large epidemics. With increased availability of sequencing it is possible for agile networks of sequencing facilities to collaborate on the sequencing and analysis of epidemic genomic data. In response to the ongoing SARS-CoV-2 pandemic in the United Kingdom, the COVID-19 Genomics UK (COG-UK) consortium was formed with the aim of rapidly sequencing SARS-CoV-2 genomes as part of a national-scale genomic surveillance strategy. The network consists of universities, academic institutes, regional sequencing centres and the four UK Public Health Agencies. We describe the development and deployment of Majora, an encompassing digital infrastructure to address the challenge of collecting and integrating both genomic sequencing data and sample-associated metadata produced across the COG-UK network. The system was designed and implemented pragmatically to stand up capacity rapidly in a pandemic caused by a novel virus. This approach has underpinned the success of COG-UK, which has rapidly become the leading contributor of SARS-CoV-2 genomes to international databases and has generated over 60,000 sequences to date.

In Silico Modeling of Virus Particle Propagation and Infectivity along the Respiratory Tract: A Case Study for SARS-COV-2 (preprint)

Respiratory viruses including Respiratory syncytial virus (RSV), influenza virus and cornaviruses such as Middle Eastern respiratory virus (MERS) and SARS-CoV-2 infect and cause serious and sometimes fatal disease in thousands of people annually. It is critical to understand virus propagation dynamics within the respiratory system because new insights will increase our understanding of virus pathogenesis and enable infection patterns to be more predictable in vivo, which will enhance targeting of vaccines and drug delivery. This study presents a computational model of virus propagation within the respiratory tract network. The model includes the generation network branch structure of the respiratory tract, biophysical and infectivity properties of the virus, as well as air flow models that aid the circulation of the virus particles. The model can also consider the impact of the immune response aim to inhibit virus replication and spread. The model was applied to the SARS-CoV-2 virus by integrating data on its life-cycle, as well as density of Angiotensin Converting Enzyme (ACE2) expressing cells along the respiratory tract network. Using physiological data associated with the respiratory rate and virus load that is inhaled, the model can improve our understanding of the concentration and spatiotemporal dynamics of virus.

SARS-CoV-2 infects brain choroid plexus and disrupts the blood-CSF-barrier (preprint)

Coronavirus disease-19 (COVID-19), caused by the SARS-CoV-2 virus, leads primarily to respiratory symptoms that can be fatal, particularly in at risk individuals. However, neurological symptoms have also been observed in patients, including headache, seizures, stroke, and fatigue. The cause of these complications is not yet known, and whether they are due to a direct infection of neural cells, such as neurons and astrocytes, or through indirect effects on supportive brain cells, is unknown. Here, we use brain organoids to examine SARS-CoV-2 neurotropism. We examine expression of the key viral receptor ACE2 in single-cell RNA sequencing (scRNA-seq) revealing that only a subset of choroid plexus cells but not neurons or neural progenitors express this entry factor. We then challenge organoids with both SARS-CoV-2 spike protein pseudovirus and live virus to demonstrate high viral tropism for choroid plexus epithelial cells but not stromal cells, and little to no infection of neurons or glia. We find that infected cells of the choroid plexus are an apolipoprotein and ACE2 expressing subset of epithelial barrier cells. Finally, we show that infection with live SARS-CoV-2 leads to barrier breakdown of the choroid plexus. These findings suggest that neurological complications may result from effects on the choroid plexus, an important barrier that normally prevents entry of immune cells and cytokines into the cerebrospinal fluid (CSF) and brain.

Covidex: an ultrafast and accurate tool for virus subtyping (preprint)

Covidex is an open-source, alignment-free machine learning subtyping tool for viral species. It is a shiny app that allows a fast and accurate classification in pre-defined clusters for SARS-CoV-2 and FMDV genome sequences. The user can also build its own classification models with the Covidex model generator. AvailabilityCovidex is open-source, cross-platform compatible, and is available under the terms of the GNU General Public License v3 (http://www.gnu.org/licenses/gpl.txt). Covidex is available via SourceForge https://sourceforge.net/projects/covidex or the web application https://cacciabue.shinyapps.io/shiny2/ Contactcacciabue.marco@inta.gob.ar; marcocacciabue@yahoo.com

Characterization of cells susceptible to SARS-COV-2 and methods for detection of neutralizing antibody by focus forming assay (preprint)

The SARS-CoV-2 outbreak and subsequent COVID-19 pandemic have highlighted the urgent need to determine what cells are susceptible to infection and for assays to detect and quantify SARS-CoV-2. Furthermore, the ongoing efforts for vaccine development have necessitated the development of rapid, high-throughput methods of quantifying infectious SARS-CoV-2, as well as the ability to screen human polyclonal sera samples for neutralizing antibodies against SARS-CoV-2. To this end, our lab has adapted focus forming assays for SARS-CoV-2 using Vero CCL-81 cells, referred to in this text as Vero WHO. Using the focus forming assay as the basis for screening cell susceptibility and to develop a focus reduction neutralization test. We have shown that this assay is a sensitive tool for determining SARS-CoV-2 neutralizing antibody titer in human, non-human primate, and mouse polyclonal sera following SARS-CoV-2 exposure. Additionally, we describe the viral growth kinetics of SARS-CoV-2 in a variety of different immortalized cell lines and demonstrate via human ACE2 and viral spike protein expression that these cell lines can support viral entry and replication.

Ubiquitous Forbidden Order in R-group classified protein sequence of SARS-CoV-2 and other viruses (preprint)

Each amino acid in a polypeptide chain has a distinctive R-group associated with it. We report here a novel method of species characterization based upon the order of these R-group classified amino acids in the linear sequence of the side chains associated with the codon triplets. In an otherwise pseudo-random sequence, we search for forbidden combinations of kth order. We applied this method to analyze the available protein sequences of various viruses including SARS-CoV-2. We found that these ubiquitous forbidden orders (UFO) are unique to each of the viruses we analyzed. This unique structure of the viruses may provide an insight into viruses chemical behavior and the folding patterns of the proteins. This finding may have a broad significance for the analysis of coding sequences of species in general.

Retrospective screening of routine respiratory samples revealed undetected community transmission and missed intervention opportunities for SARS-CoV-2 in the United Kingdom (preprint)

In the early phases of the SARS coronavirus type 2 (SARS-CoV-2) pandemic, testing focused on individuals fitting a strict case definition involving a limited set of symptoms together with an identified epidemiological risk, such as contact with an infected individual or travel to a high-risk area. To assess whether this impaired our ability to detect and control early introductions of the virus into the UK, we PCR-tested archival specimens collected on admission to a large UK teaching hospital who retrospectively were identified as having a clinical presentation compatible with COVID-19. In addition, we screened available archival specimens submitted for respiratory virus diagnosis, and dating back to early January 2020, for the presence of SARS-CoV-2 RNA. Our data provides evidence for widespread community circulation of SARS-CoV2 in early February 2020 and into March that was undetected at the time due to restrictive case definitions informing testing policy. Genome sequence data showed that many of these early cases were infected with a distinct lineage of the virus. Sequences obtained from the first officially recorded case in Nottinghamshire - a traveller returning from Daegu, South Korea - also clustered with these early UK sequences suggesting acquisition of the virus occurred in the UK and not Daegu. Analysis of a larger sample of sequences obtained in the Nottinghamshire area revealed multiple viral introductions, mainly in late February and through March. These data highlight the importance of timely and extensive community testing to prevent future widespread transmission of the virus.

Targeting pentose phosphate pathway for SARS-CoV-2 therapy (preprint)

It becomes more and more obvious that deregulation of host metabolism play an important role in SARS-CoV-2 pathogenesis with implication for increased risk of severe course of COVID-19. Furthermore, it is expected that COVID-19 patients recovered from severe disease may experience long-term metabolic disorders. Thereby understanding the consequences of SARS-CoV-2 infection on host metabolism can facilitate efforts for effective treatment option. We have previously shown that SARS-CoV-2-infected cells undergo a shift towards glycolysis and that 2-deoxy-D-glucose (2DG) inhibits SARS-CoV-2 replication. Here, we show that also pentose phosphate pathway (PPP) is remarkably deregulated. Since PPP supplies ribonucleotides for SARS-CoV-2 replication, this could represent an attractive target for an intervention. On that account, we employed the transketolase inhibitor benfooxythiamine and showed dose-dependent inhibition of SARS-CoV-2 in non-toxic concentrations. Importantly, the antiviral efficacy of benfooxythiamine was further increased in combination with 2DG.

Temporal landscape of mutation accumulation in SARS-CoV-2 genomes from Bangladesh: possible implications from the ongoing outbreak in Bangladesh (preprint)

Along with intrinsic evolution, adaptation to selective pressure in new environments might have resulted in the circulatory SARS-CoV-2 strains in response to the geoenvironmental conditions of a country and the demographic profile of its population. Thus the analysis of genomic mutations of these circulatory strains may give an insight into the molecular basis of SARS-CoV-2 pathogenesis and evolution favoring the development of effective treatment and containment strategies. With this target, the current study traced the evolutionary route and mutational frequency of 198 Bangladesh originated SARS-CoV-2 genomic sequences available in the GISAID platform over a period of 13 weeks as of 14 July 2020. The analyses were performed using MEGA 7, Swiss Model Repository, Virus Pathogen Resource and Jalview visualization. Our analysis identified that majority of the circulating strains in the country belong to B and/or L type among cluster A to Z and strikingly differ from both the reference genome and the first sequenced genome from Bangladesh. Mutations in Nonspecific protein 2 (NSP2), NSP3, RNA dependent RNA polymerase (RdRp), Helicase, Spike, ORF3a, and Nucleocapsid (N) protein were common in the circulating strains with varying degrees and the most unique mutations(UM) were found in NSP3 (UM-18). But no or limited changes were observed in NSP9, NSP11, E (Envelope), NSP7a, ORF 6, and ORF 7b suggesting the possible conserved functions of those proteins in SARS-CoV-2 propagation. However, along with D614G mutation, more than 20 different mutations in the Spike protein were detected basically in the S2 domain. Besides, mutations in SR-rich region of N protein and P323L in RDRP were also present. However, the mutation accumulation showed an association with sex and age of the COVID-19 positive cases. So, identification of these mutational accumulation patterns may greatly facilitate drug/ vaccine development deciphering the age and the sex dependent differential susceptibility to COVID-19.

Evaluating the effects of SARS-CoV-2 Spike mutation D614G on transmissibility and pathogenicity (preprint)

Global dispersal and increasing frequency of the SARS-CoV-2 Spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of Spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large data set, well represented by both Spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the Spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant.